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Abstract 



We use the Ricci flow with surgery to study four-dimensional SU{2) x C/(l)-symmetric metrics 
on a manifold with fixed boundaiy given by a squashed 3-sphere. Depending on the initial metric 
we show that the flow converges to either the Taub-Bolt or the Taub-NUT metric, the latter case 
J> ' potentially requiring surgery at some point in the evolution. The Ricci flow allows us to explore 

^ . the EucUdean action landscape within this symmetry class. This work extends the recent work of 

\^ I Headrick and Wiseman [Sj to more interesting topologies. 

o 1 Introduction 

i> ■ 

^ ' Ricci flow is a geometrie flow within the space of Riemannian metrics on a given n-dimensional 

manifold Ai defined by: 

^1 Tt9it,x) =-2RiCg^t,x) ^^^^ 

g{0,x) =go{x) 

where g{t) is the metric at flow time t and RiCg(^t) is its associated Ricci tensor. On a closed manifold, 
the PDE dl.ll ) admits short time existence 161, but will generally exhibit finite time blow-up. Work 
of Hamilton and Perelman has provided a means of continuing the flow beyond such a singularity, 
namely Ricci flow with surgery. There has recently been much interest in the Ricci flow with surgery, 
especially since the publication by Perelman of a proof of the Poincaré conjecture based on this 
method ifTOlfTTI. 
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For many practical purposes it is more convenient to consider a slightly modified Ricci flow, 
defìned by: 

Ìg{t,x) =-2i?ic,(,,,) + 2(Ve(t))s ^^^^ 
g{0,x) =go{x) 

Here ^{t) is an arbitrary vector field on A4 and indicates that the tensor enclosed should be 
symmetrized on ali its indices. The flows (11.11) and (ll.2b are equivalent if we consider the flow to 
act on the space of metrics modulo dijfeomorphisms as the second temi in (11.21 ) on its own would 
define a flow through the equivalence class of a metric under diffeomoiphisms. The advantage of 
(11.21 ) however is that one can pick the gauge (t), i.e. a diffeomorphism at each point of time, such 
that the equations become strongly parabolic. Short-time existence then foUows from standard results 
|[T]|2l. Moreover the stabili ty of numerical codes is improved in this formulation of the Ricci flow. 

In this paper we are going to study the Ricci flow on a class of four-dimensional manifolds Ai 
whose boundary is given by a squashed 3-sphere. The physical motivation to study infiUing metrics 
for a manifold with fixed boundary metric derives from an attempt to understand the thermodynamics 
of general relativity. In the Euclidean approach to the problem the metric is analytically continued to 
the Riemannian sector and - in analogy with statistical mechanics - a partition function associated 
with the canonical ensemble is defined by 

Z{g)= [ d[<7]e-^[^]. (1.3) 
Jm 

Here S [g] is the Euclidean action (12.41) and the integration is performed over ali metrics on Ai ad- 
mitting the same fixed boundary metric but (perhaps topologically) different interiors. The fixed 
boundary metric can be understood in a thermodynamic sense as holding the system at fixed temper- 
ature. However, the integrai (11.31 ) is not well-defined in general and in the one-loop semi-classical 
approximation renormalization schemes bave to be employed to extract the correct values for the en- 
tropy and other thermodynamic quantities from the partition function ||5l. In practical evaluations of 
(11.31 ) one usually restricts to certain symmetry classes and invokes a saddle point approximation: the 
main contributions to the integrai (11.31 ) will come from extremal points of the action, which are found 
to be Ricci-flat metrics. In this setting it is an interesting mathematical problem on its own to classify 
the number of infilling Ricci-flat metrics for a given boundary metric. This question has so far been 
answered analytically only in restricted symmetry classes (like the one studied in this paper). On the 
other band, being a gradient flow, the Ricci flow provides a useful tool to explore the action landscape 
of Riemannian metrics. For metrics less constrained by symmetry assumptions, the Ricci flow may 
be thought of as a relaxation technique for finding new Ricci-flat metrics. 

We direct the reader interested in further details, especially on the relation of the Ricci-flow to the 
thermodynamics of black holes and the renormalization group in quantum field theory, to [i8J. In that 
paper the authors study the canonical ensemble for gravity in a box, modelled by four-dimensional 
manifolds with a boundary of topology x S^Q Furthermore, they restrict to metrics admitting a 
U{1) X SO (3) action by isometryH At high temperatures, i.e. a small enough size of the direction, 
they find three infilling Ricci-flat geometries, which correspond to saddle points of the action: two 

'The direction arises from the periodic time coordinate of the Lorentzian metric after analytic continuation. 
^Note that symmetry in the initial data is preserved under the Ricci flow. 
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Schwarzschild black holes of different horizon areas (with topology B'^ x S"^) and hot fìat space 
(topology X B^). Considering perturbations of the small Schwai^zschild solution it is fìnally shown 
numerically how the Ricci-flow converges to either the large black hole or hot fiat space, with the 
latter case enforcing a surgery at some point. 

The results of the present paper provide a naturai extension of the work of HI to other topolo- 
gies. Here the boundary will be given by a squashed and we will restrict attention to manifolds 
endowed with metrics admitting an SU{2) x U{1) symmetry, so called biaxial Bianchi IX metrics. 
The Ricci fiat metrics, i.e. the fìxed points of the flow (11.11) . within this class are the one parameter 
families of Taub-Bolt and Taub-NUT metrics, which differ in their topology. Taub-Bolt is a metric on 
CP^ \ {open ball}, whereas Taub-NUT is defìned on B'^. We will discuss some properties of these 
metrics, which play a role in Euclidean quantum gravity Q, in section 2. We also establish that fixing 
the squashing parameter of the 3-sphere on the boundary in a certain range in fact singles out three 
infilling Ricci fiat geometries from these families: Two Bolt- and one NUT- solution. We demon- 
strate, by finding a negative eigenvalue of the linearized operator, that one of these Bolt solutions is 
unstable under the Ricci-flow. The full non-linear evolution of the unstable solution is finally studied 
numerically. We find that the unstable Bolt metric will flow to either the second (stable) Bolt metric 
or the NUT metric, depending on how the metric is perturbed initially. Geometrically the perturbation 
conesponds to whether one slightly increases or decreases the ai^ea of the small Bolt minimal surface. 
The cruciai and most interesting feature, however, is that the latter flow from the Bolt to the NUT 
solution requires surgery because the NUT and Bolt solutions bave different topology, as mentioned 
above. Hence we will explain our method of controlling the curvature blow up and how the surgery 
is finally performed. The behaviour under perturbations just described might be expected from an 
analysis of the gravitational Euclidean action for the solutions as illustrated below. 

This work, together with that of Headrick and Wiseman, provide concrete examples of the two 
types of surgery required by the Ricci flow in 4 dimensions as anticipated by Hamilton in fTI, which 
locaUy look like 

S^xB^^B^yiS^ and S"^ x B"^ ^ B^ x . (1.4) 

respectively. These surgeries ai^e expected to be iiTcversible. 

To our knowledge the question of existence and uniqueness for the Ricci flow on manifolds with 
fixed boundary metric has not been answered in any generality. The reason is that in the strongly 
parabolic formulation (11.21) one typically ends up with mixed Dirichlet and Neumann conditions at the 
boundary. We will therefore prove uniqueness of the evolution for our Ricci-flow equations directly 
in appendix A. 

2 NUTs and Bolts 

2.1 Regular Bìaxìal Bianchì IX Metrics 

The biaxial Bianchi-IX metrics are a class of metrics admitting an SU (2) lxU {1)r symmetry group. 
They may be written in terms of the usuai left-invariant su{2) one-forms cjj in the form: 

ds"^ = airfdr"^ + b{rf {al + al) + c(r) V| (2. 1) 

The SU{2)l symmetry arises from the left action of SU {2) under which the one-forms ai are by 
definition invariant. In addition, since ai and a2 appear in the metric with the same coefficient there 
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is a right action by a U{1) subgroup of SU {2) which acts by isometries. We will be interested in 
manifolds with a boundaiy, so the radiai coordinate r will take values in the range < r < ri. The 
boundary at ri is topologically a 3-sphere and the induced metric that of a homogeneously squashed 
S^. We assume that the functions a, h, c are smooth and non zero on the interval and further that 
{r = vq} is a fixed point set of the U{1)r action. This means that c(ro) = and there will in general 
be a singularity at r = unless certain conditions are met. Without loss of generality, we assume 
that ro = and a(r)^ ~ Oq + O(r^). If the metric is to be regular at r = there are two possibili ties 
Bl: 

Case 1 6(0) 7^ 0. In this case regularity requires that 

<rf ~ + Oir^) (2.2) 

as r — ^ 0. A coordinate transformation then shows that r = tq is a coordinate singularity, similar to 
the origin of in polar coordinates. The point set r = is then a metric 2-sphere known as a Bolt. 
Geometrically it is a minimal surface rather than a boundary. 
Case 2 6(0) = 0. In this case regularity requires that 

b{rf ~ c{rf ~ + O(r^) (2.3) 

as r — > tq. a coordinate transformation shows that r = is a coordinate singularity, similar to the 
origin of in spherical polar coordinates. The point set r = is then a point known as a Nat. 



2.2 Ricci Fiat Metrics 

Clearly from (ll.lb the fìxed points of the Ricci flow are metrics satisfying Ricg = 0. Such metrics are 
sometimes known as gravitational instantons and occur as extremal points of the Euclidean Einstein- 
Hilbert action: 

^[5] = -TTT^ I dV^R--^ [ dSVhK. (2.4) 

With a carefully chosen metric on the space of Riemannian metrics, the Ricci flow may be thought of 
as a gradient flow of this action, as discussed by Headrick and Wiseman [SJ. 



2.2.1 Taub-Bolt 

In order to find Ricci flat infilling metrics for a squashed boundary, we consider Ricci flat biaxial 
Bianchi IX metrics whose radiai slices are squashed S^. There are two families of such metrics 
con^esponding to the two possible regular structures at the fixed points of the U{1) action. The Taub- 
Bolt family of metrics is given by: 

ds' = ^4^^rfr^ + {r' - n') {ai + aj) + ('-^^P^] al (2.5) 

This metric is Ricci flat and regular for r in the range 2n < r < 00 and there is a bolt at r = 2n. 
In order that this metric may be considered as infilling a squashed we introduce a boundary by 
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restricting r to 2n < r < R. It will be convenient to change the radiai coordinate and introduce new 



parameters a = y2 + ^,P = \/3n instead of n, R so that the metric instead takes the form: 

j 2 oq2 2 ~ P ì 1 2 , 02 \^ ~ P ì { 2 , 2\ , <^n2 2 P 2 ^\ 

= 8/5 « (^2_2p2)4 ^^ + /5 (^2 _ 2p2)2 (^1 + ^2) + 2/3 a ^tT^^s (2-6) 



Now the radiai coordinate is in the range < p < 1. For > 2 this is a regular metric on 
a manifold with boundary which is topologically CP^ with an open 4-ball removed. In the limit 
^ 2 the boundary receeds to infinity and we bave the asymptotically locally fìat Taub-Bolt metric 
on - {pt.}. 

We seek Ricci fìat metrics which bave a given metric on the boundary S^, taken to be: 

h = p^{al + al)+u^al (2.7) 



h can only be the induced metric on the boundary of a Taub-Bolt solution for certain values of = 
I [x^, which we refer to as the squashing pai^ameter. Ifr^ > 1 — |33 + §33 =t^ then one cannot 
find a Taub-Bolt solution whose boundary has metric (12.71 ). For < there are two choices of 
(q, /9) such that the boundary metric is h. These can be distinguished via the parameter /9 which 
parameterises the size of the bolt. The two solutions ai^e thus distinguished by having a 'small bolt' 
or 'large bolt'. In the criticai case where r = Tc, the small and lai^ge bolt solutions coincide and we 
have only one infìUing Taub-Bolt geometry. 



2.2.2 Taub-Nut 

We now turn to the other family of Ricci fiat biaxial Bianchi IX metrics which may be considered as 
infìlling geometries for the boundary 5^ with metric (12.71 ). These contain a nut and are known as the 
Euclidean self-dual Taub-NUT solutions. In their usuai form they have metric: 

^^2 ^ l±J^^^2 4„2!:^^2 / 2 _ ^2^ / 2 2N _ (2.8) 

This metric is regular and Ricci fiat for r in the range n < r < oo. At r = n the metric has a regular 
nut and we introduce a finite boundary by restricting r to n < r < R. As in the previous case, it will 
be convenient to change to a different radiai coordinate and new parameters 7 = -^\/2n, S = \/R so 
that the metric takes the form: 

ds^ = (p2 + ^2)^^2 ^ ^2^2 j^(^2 ^ ^2) (^2 ^ ^2) ^ -^^af) , (2.9) 

where the radiai coordinate is in the range < p < 1. One can see by inspection that this metric 
satisfies the conditions for p = to be a regular nut. This gives a regular metric on B^, the closed 
4-ball. There is exactly one Taub-Nut metric which infiUs the squashed (12.71 ) for r in the range 
< r < 1. 



2.2.3 Action 

As we know that the Ricci flow is a gradient flow of the Euclidean action (12.41 ). it is convenient to plot 
the action of the Ricci fiat solutions. We fix the volume of the boundary metric and plot the difference 
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Figure 1 : A plot of 5'Boit — 'S'Nut as a function of boundary squashing 

SboM — ^Nut against r for both branches of the infìlling Taub-Bolt solutions. This plot is included 
as Figure [T] We note that the solution with the greatest action is always the small bolt solution. The 
solution with the least action is either the large bolt or the nut, according a.s t < thp or t > thp- At 
T = thp there is a Hawking-Page type transition as the global minimizer of the action changes from 
one solution to the other 

In many ways, the Taub-Bolt geometries are similar to the Euclidean Schwarzschild geometries 
considered in IH. After analytic continuation, the horizon of the Schwarzschild black hole becomes 
a regular bolt provided the Wick rotated time coordinate has some suitable period. If we seek ge- 
ometries which infili an 5^ x boundary, then we fìnd that this can only be done with a Euclidean 
Schwarzschild metric for a particular range of the ratio 13/R. As with Taub-Bolt there is a criticai 
value for this ratio where there is a bifurcation and two solutions appear. The Schwarzschild solution 
with a small bolt has greatest action and is unstable under the Ricci flow, flowing to either the large 
bolt solution or to periodically identified fiat space. 

Motivated by this analogy to Euclidean Schwarzschild, we expect the small bolt solution to be 
unstable under Ricci flow in the case of biaxial Bianchi IX symmetry. More precisely, perturbing the 
small bolt solution the metric should flow to either the large bolt solution or the Taub-Nut solution, 
the latter flow requiring a topology changing surgery. We will show in the remainder of the paper that 
the flow indeed exhibits the behaviour described. 
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3 Ricci flow simulations 



3.1 Bolt topology - the linearized problem 

We will begin our simulations by perturbing the small bolt solution as we expect this to be unstable 
against the Ricci flow. For the flow on a manifold with bolt topology we make the metric ansatz: 

ds^ = e2^(^'*)dr2 + e2^('^'*)(<T? + aj) + r'e^''^^^'^ al (3.1) 

with < r < 1. The Taub-Bolt solutions may be cast into this form. We also make a choice of ^, the 
vector field which specifies the gauge in (11.21 ). We take the usuai {9, (p, ip) coordinates for SU (2) and 
set 

^ = [-{g^„ irdi. dx") + (g^,^ * d * dx'')o] dx^" = -{2B' + C - A')dr, (3.2) 

where the subscript indicates that the term should be evaluated on some fixed reference metric, 
which we take to be (12.61 ). One can check that ^ vanishes identically in these coordinates. These 
are the coordinates singled out by the gauge choice, which breaks the diffeomorphism invariance of 
the Ricci flow equations. This choice of gauge, similar to that chosen for the 'deTurck trick' ||2l 
makes the Ricci flow equations strongly pai^abolic. We wish to evolve the functions on the interval 
< r < 1, however r = is not a boundary for the manifold, but a point at which the coordinates 
break down. This is reflected in the equations we derive from (13.11 13.21 ) as some coefficients of the 
PDE are not bounded as r ^ 0, cf. the system (IA.2I )- (IA.4I ) in the appendix. The PDE is expected 
to admit a well posed initial value formulation, provided that we impose that the functions A, B,C 
extend smoothly to even functions on the interval — 1 < r < 1 as will be discussed in appendix A, 
where we indeed prove the uniqueness of solutions with given initial and boundary data. Doubling 
the interval is therefore a convenient strategy to remove the apparent boundary at r = 0. The pseudo- 
spectral collocation methods yj we bave used for the numerical simulations benefit to a large extent 
from this idea. We now need three boundary conditions at r = ijl Two of these arise by fixing the 
metric on the boundary. For well posedness we require a third, which comes from requiring the vector 
field ^ to vanish at the boundary. This ensures that the boundary remains at r = 1. The boundary 
conditions are then: 

B{l,t) = 5(1,0) 

C{l,t) = C{1,0) (3.3) 
e(l,t) = 

We expect that the small bolt solution will be unstable under this flow and that the large bolt 
solution will be stable. In order to check this, we linearize the system of equations obtained above 
about the known Taub-Bolt solution and seek solutions of the form 

A{r,t)=Ao{r)+e-^'air), B{r,t) = Bo{r) + e~^%{r), C{r,t) = Co{r) + e'^'cir), (3.4) 

where a, b, c are considered small. It can be shown that A is then an eigenvalue of the Lichnerowicz 
Laplacian. If the Lichnerowicz Laplacian has a negative eigenmode, the perturbations will grow 

■'those for r = — 1 follow by symmetry 
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Figure 2: A plot of A^m as a function of a 

without bound and the solution will be unstable. We bave calculated numerically, using a pseudo- 
spectral metbod based on collocation in Cbebysbev points llT3l Cbap. 9] tbe lowest eigenvalue of tbe 
Licbnerowicz Laplacian, in tbe gauge ^. It is plotted in Figure (|2]l as a function of a, for 0^ = 3. 
We tìnd, as expected, tbat tbe lowest eigenvalue is negative for < wbicb is tbe range of 
parameters corresponding to tbe small bolt. Tbe eigenvalue vanisbes at tbe criticai point wbere tbe 
small and large bolt coincide and becomes positive for tbe large bolt. Witb 38 sampling points, tbe 
cbange in sign occurs at tbe tbeoretical value of al to witbin 10"^^. Furtbermore, as a ^ 2, tbe 
boundary approacbes infinity and we recover tbe result of Young tlSil tbat ALE Taub-Bolt bas a 
negative Licbnerowicz mode witb A G (—0.200, —0.201)0 Tbe cboice of [ì above is so as to make 
a direct comparison witb Young 's work possible. See also |[T4l wbere a cosmological Constant is 
included. Sucb negative modes may be interpreted as indicating a semi-classical instability of tbe 
metric considered as a gravitational instanton |51. Tbey also indicate a dynamical instability witbin 
general relativi ty of tbe 5-dimensional ultra-static vacuum metric obtained by adding a — dr^ term to 
tbe metric. 

^For values of near 2, one must increase the number of sampling points to obtain accurate results. 
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Figure 3: Plots of Vi{x) ^nd V2(x) at selected times during the flow for (3 = l,a^ = 2.5. The flow 
proceeds in the direction of the arrow. The two infilling bolt geometries are shown as solid lines. 



3.2 Taub-Bolt Topology - the Initial Flow 

We use the solutions to the linearised problem to seed the flow for the full non-linear equations. We 
consider the flow with initial conditions: 

A{r) = AQ{r) +ea-{r), B{r) = Bo{r) + eb^{r), C(r) = Co(r) + ec_(r) (3.5) 

Where the functions a_ etc. are the negative mode found for the linearised problem, normalised so 
that 6_(0) = 1. In this case a perturbation with e > increases the area of the bolt, while one with 
e < decreases it. 

We first consider perturbing the small bolt solution with e > which produces a small increase 
in the area of the bolt. The PDEs were again discretized in space using pseudo-spectral collocation in 
Chebyshev points. The time-stepping was performed with MATLAB's odel 5s numerical integrator. 
We find that for a range of values of and e > the metric at time t approaches the large bolt solution 
as t becomes large. In order to plot the results, we first cast the metric in the gauge invariant form: 

= dx^ + Viixf {af + al) + xV2(x)Vi (3.6) 

and then plot the functions Vi and V2. We see that the flow starting at the small bolt converges to 
the large bolt. Taking 102 spatial sampling points and (3 = l,a^ = 2.5, we find that for t > 40 the 
functions bave uniformly approached the large bolt solution to within 10~^, with better accuracy for 
more sampling pointsjf] For larger values of we also get better accuracy, indeed for = 4 we can 
achieve 10^^^ accuracy with only 42 sampling points. 

We now turn to the situation where e < 0. In this case, the numerical simulations appear to exhibit 
a finite time blow-up which occurs because the bolt becomes smaller and smaller and approaches zero 
size in finite time. This behaviour where an 5^ collapses to a point in finite time is chaiacteristic of 
the Ricci flow. In order to continue the flow, we must perform a surgery, turning the bolt at the origin 
to a nut. To decide when such a surgery should he performed, we monitor the value of the coordinate 
invariant scalar K = Riem'^\r=o- We stop the flow when t = t^, defined by K{tc) > 10^. 

^However, there is a limit due to the fact that for very large N the matrices typically become ill-conditioned. 
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Figure 4: Plots of Vi (x) at various times, both before surgery (on the left) and after surgery (on the right). 
The solid red line shows the function after surgery on both plots. The solid black lines are the functions 
for Taub-Bolt on the right and Taub-NUT on the left. 



3.3 Taub-NUT topology - Surgery and Beyond 

The idea behind surgery for the Ricci flow is that given a solution to the Ricci flow equations ex- 
hibiting finite time blow-up at a time T, we should take the Riemannian manifold at a time T — 6 
and remove a small region of size e from the manifold and glue in a new piece which has a different 
topology and smoothly extends the metric on the rest of the manifold. We then re-stait the Ricci flow. 
If this is done appropriately, the resulting flow should not depend on the details of the gluing in the 
Ihnit 6^0. 

In our case, the minimal 5^ at r = appears to be shrinking to zero size in a finite time. In order 
to continue the flow beyond this singulaiity, we remove the region r < e from the manifold and glue 
in a 4-ball. This means that the new metric will take the form: 



ds' = e'^ ^dr' + [e'^'^ial + ai) + ai) j (3.7) 

We require that this agrees with the metric at the point we stopped the flow, g{tc), for r > e and is 
every where smooth. This leaves considerable ambiguity in prescribing the metric functions for r < e, 
however numerica! experimentation showed that the long term properties of the flow ai^e insensitive 
to the precise gluing technique. 

We now wish to continue the Ricci flow with the new initial conditions given by the metric func- 
tions after surgery. We will again introduce a gauge field ^ to make the Ricci flow equations strongly 
parabolic. We make a similar choice of gauge, given by: 

^= [-{g^^i.di.dx'') + {g^^i.di.dx'')o]dx^' = - \2À' + 2{r'By-3{r'C)'] dr. (3.8) 

This time we use the Taub-NUT metric in coordinates (12.91 ) as a reference metric. The equations 
which arise are once again symmetric under r —r and so we can set the boundary conditions at the 
origin by requiring that the functions A, B,C extend smoothly to an even function on [—1, 1]. The 
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Figure 5: Plots of xV2(x) at various times, both before surgery (on the left) and after surgery (on the 
right). The solid red line shows the function after surgery on both plots. The solid black lines are the 
functions for Taub-Bolt on the right and Taub-NUT on the left. 

symmetric boundary conditions are then given by: 



which indicate that the boundary metric is fixed and the diffeomorphism ^ fixes r = 1. This final 
condition is not satisfied by the functions obtained via surgery, however a near identity change of the 
r coordinate can make the initial conditions compatible with the boundary conditions. 

We can show (see Appendix) that the system of PDEs obtained from (11.21) with the gauge choice 
(13.81) and boundary conditions (13.91 ) has the uniqueness property. We believe that short temi existence 
also holds, making the problem well posed. One might think that the parity condition could be 
replaced by ^'(0, t) = etc. however it appears that with this boundary condition the system is 
ili posed and short term solutions do not exist for generic initial data. 

The Ricci flow equations were solved numerically both before and after surgery using the same 
technique as for the flow not requiring surgery, starting from a Taub-Bolt metric with /3 = 1, = 3.5. 
The results are shown in figures|4]and[5] For these plots it is more convenient to plot Vi and x^2(x)' 
where the functions Vi are defined as in the previous section. We see that although the curvature blows 
up at r = under the flow within the bolt topology, the gauge invariant functions are well behaved 
and the change due to the surgery is slight. We find that for sufficiently many sampling points the 
final functions approach the functions of Taub-NUT uniformly to within 10"^. Experimentation with 
arbitrary initial conditions suggests that the Taub-NUT metric is the final state of the flow for any 
metric within the B'^ topology with appropriately squashed S"^ boundary conditions. 



{A-C + B){l,t) 

(I-CKi,t) 







(A-C + S)(1,0) 
(I-C)(1,0) 



(3.9) 
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4 Conclusion and Open Questions 



In this paper we explored the Euclidean action landscape of biaxial Bianchi IX metrics using the Ricci 
flow with surgery. We demonstrated the existence of three saddle points of the Euclidean action and 
showed numerically how the Ricci-flow takes us from one to the other, namely from the unstable 
Bolt-solution to either the stable Bolt or the NUT-solution depending on the perturbation. 

It is an ambitious analytical problem to make the observations of the present paper rigorous, which 
would involve applying the methods of Hamilton and Perelman to study the long-time behaviour of 
the Ricci-flow. We finish the paper with two conjectures, which we hope to address rigorously in the 
future. The first expresses our belief that we bave exhausted the quaUtative features of the Ricci-flow 
within this symmetry class by our study: 

Conjecture 4.1. For arbitrary biaxial Bianchi-IX symmetric initial-data and fixed metric on the 
boundary given by Ii2. 7ÌI with < r < 1, Ricci-flow plus at most one surgery will converge to 
a member of either the Taub-Bolt or the Taub-NUT family. 

Restricting to the subcase of NUT-topology, one might - encouraged by the final remark of the 
previous section - hope to prove the stronger 



Conjecture 4.2. For any biaxial Bianchi-IX symmetric initial-data with Taub-NUT topology ( 



the Ricci-flow admits long-time existence and converges to a member of the Taub-NUT family. 

It should be interesting to see how the different geometrie features of the Bolt and the Nut make 
their way into the estimates necessary to establish the long-time behaviour of the flow. 

Let US finally mention possible extensions of the present work. With the tools at band it should 
be possible to include a cosmologica! Constant into the formulae and impose AdS boundary condi- 
tions. We expect similai^ results to hold with Taub-Bolt and Taub-NUT being replaced by their AdS 
generalizations. 

Moreover, there is a topological generalization if one considers the more general Lens spaces 
instead of the squashed to function as the fixed boundary metric. For the case of MP'^ we expect 
the Eguchi-Hanson metric to be a stable infilUng solution. 
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A Uniqueness 

In this appendix we are going to prove uniqueness of the solution to our modified Ricci-flow equa- 
tions. As we noted in the introduction, short-time existence and uniqueness has been proven for the 
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Ricci flow in general on closed manifolds. On the other band, for manifolds with boundary the only 
general existence and uniqueness result we know of is [12J where the Neumann problem 



dt9ij {x, t) = -2Rij {x, t) 
9ij ix,0) = gij{x) (A.l) 
Kah = ondW 

with Kah being the extrinsic curvature on the boundaiy, is studied. It tums out that the boundary has 
to satisfy certain geometrie conditions in order for (lA.ll ) to admit short-time existence: For a totally 
geodesie boundary existence and uniqueness can be shown |[T2l by an application of the general results 
on parabolic boundary value problems described in IH. 

The Dirichlet problem, the problem of fixing the metric itself on the boundary, has not yet been 
considered in generality. The reason seems to be that the deTurck technique will always introduce a 
derivative condition at the boundary and the resulting boundary conditions for the modified flow are 
of non-standai^d type. Note that our boundary condition (I3.3b is typical in that sense. 

This motivates a more detailed study of existence and uniqueness within the biaxial symmetry 
class considered in the paper. Here we will prove uniqueness of the modified Ricci-flow. A more 
practical reason which originally initiated the study was the occurrence of numerical solutions blow- 
ing up near the origin in our evolution. It turns out that in addition to the boundary conditions a 
topological condition has to be imposed near the origin con^esponding to the avoidance of a conical 
singularity. Without this additional condition our uniqueness proof would not go throughH 
The Ricci-flow equations for the metric ansatz (13.11 ) constitute the following non-linear system 

A = e-2A^" _ Ìe-2A^' + ^e"2^C" + + 2i?'2 - ^'2) (A.2) 

B = • + -e-2^S' - + (A.3) 

r 2 

C = e-2^-C" + -e-2^C"--rV^^+2^ (A.4) 
r 2 

The vaiiable r ranges from to 1. We will prove uniqueness within the space of smooth even func- 
tìons. In particular, the Neumann conditions, A'{t, 0) = 0, B'{t, 0) = 0, C'{t, 0) = hold at r = 0. 
At r = 1 the Dirichlet conditions, B{t, 1) = Bq{1), C{t, 1) = Co(l), and the Neumann condition, 
{2B' + C — A') {t, 1) = is imposed, as was outlined in section 2. Moreover, to avoid a conical 
singularity at the origin, we must bave lim^^o — C) = log 2. 

To avoid the coupled boundary condition we make the transformation 

E = 2B + C - A (A.5) 



^Uniqueness is likely to be violated if that condition is not imposed, as illustrated by the following lower dimensionai 
elliptic example. The infilling Ricci-fìat geometries for a manifold with fìxed metric on a 5^ boundary are the fìat disc but 
also the cone. This type of non-uniqueness is expected to carry over to the parabolic case. 
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to find the system 



^4B+2C-2Eq 
^4B+2C-2E ^ 



B" + -B' - e^B+2C-2E _^ \^2^iC-2E 
r 2 

r 2 
E" + -E' +(--E' + -B' 



(A.6) 
(A.7) 
(A.8) 



C'2 + - [2B' + C'- Ef) + ir^e^^-^^ - 26^^+2^-2^ 



where now £" = at the boundary r = 1. Next we assume the existence of two even solutions 
{B, C, E) and (^B, C, E^ , both bounded with ali derivatives up to order k, say, in [0, 1] x [0, T] for 
fixed T. We are going to consider the equations for their difference 



The equations are 

^4B+2C-2Ej^ 
AB+2C-2Ep 



AB+2C-2E 



s 



iB,C,£) = [B - É,C -C,E - É 



(A.9) 



^// _|_ _|_ ^ /4B+2C-2E _ ^AB+2C-2E 



^2B+2C-2E _ ^2B+2C-2EJ _|_ ^^2 ^iC-2E _ ^AC-2E 
C" + \c' + è [e^B+2C-2È _ g4B+2C-2i?^ 
_Ì^2 ^fAC~2E _ ^AC-2È 



(A. 10) 



(A. 11) 



£" + +[--£' + -B' \+E{ e4S+2C^-2£^ _ e4B+2C-2E 



c + c 



e + 2 



B + B 



B' 



with the initial data 



2B + C-E + 2B + C-E 



[2B' + C' - £') 



1 



_| ^2 ^^iC-2E _ ^4C-2E^ _ 2 ^g2_B+2C-2£; _ ^2B+2C-2E 



5(0, r) = C(0,r) = £{0,r) = 

and boundary conditions 

B{t,l) =C{t,l) = and 1) = 

as well as the conical condition 

- 25 ~ at the origin 



(A.12) 

(A. 13) 

(A. 14) 
(A. 15) 
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plus regularity conditions at the origin arising from the fact that the functions are even. 
We now multiply equation (lA.lOI ) by (^rB + rB^ , (lA.llI ) by rC, (IA.12I) by {£ - 2B), add the terms 
up and integrate over time and space. Here fc is a fìxed Constant which we are going to determine 
below (k = 16 will do for instance). On the left band side we obtain for the B equation dA.lOI ) when 

multiplied by [rB + rB] and integrated 



t2 /-l 



ti Jo 



^AB+2C-2E QB + B 



drdt 



+ 



1 

re 
Ì2 /■! 







t=t2 
t=ti 



ti JO 



re^B+2C^2E^2^^^^ 



1 

+ 2 



t2 /-l 



il JO 



re 



AB+2C-2E -2C + 2È] B^drdt 



(A. 16) 



for the C equation (lA.l 11 ) when multiplied by Cr and integrated we find 

ft2 /■! 



/ ' / r [e4S+2^-2^ (CdtC)] drdt 
Jti Jo 



re 



4B+2C-2Ep2 



C'dr 



+ - 



t2 /■! 



2 Jti Jo 



re 



AB+2C-2E 



-AB -2C + 2E]C 



t=t2 
t=ti 

drdt 



(A. 17) 



and finally for the left band side of the £ equation ( IA.12I ) after multiplication by {£ — 2B) and 
integration 

rt2 /■! 



,'iB+2C-2E 



(dtS) {£ - 2B)] drdt 



2k 



re 



AB+2C-2E (c2 



- A£B) dr 



+ 



t2 ri 



ti JO 



2k 



-re 



4B+2C-2E 



-4È -2C + 2È) 2 _ 4f 



+ 



ti Jo 



1 1 



AB+2C-2E 



2£B] 



t=t2 
t=ti 

drdt 
drdt 



(A. 18) 



The spacetime integrals of equations (IA.16I ). (IA.17I ). (IA.18I ) not containing a time-derivative of B, C 
or £ aie controlied by 



K 



t2 /■! 



ti JO 



re 



AB+2C~2E (n2 



{B^+C^+£^) drdt 



(A. 19) 



where K is a large Constant. This follows from the fact that we control the time-derivatives of the 
functions A, B and C by some Constant. Note also that the term in the last line of (lA.lSb can be 
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bounded by 



t2 /-l 



il JO 



1 



drdt < 



t2 /■! i 
ti JO >^ 



4B+2C-2E 



8^ +B' 



drdt (A.20) 



Let US turn to the terms which are produced on the right band side when equation (lA.lOI) is multipHed 
by [rB + ré] foUowed by a spacetime integration. For the highest order derivative terms we obtain 



i2 /■! 



il JO 



B" + -B' \ (B + B] rdrdt 
r 



Ì2 



(B') rdr 



r(B'B + B'B 

rt2 /■! 



dt 



r=0 



ti JO 



r (B') drdt 



(A.21) 



The boundary terni vanishes because B satisfies Dirichlet conditions at r = 1 and Neumann conditions 
at r = (even function). Analogously we estimate the highest order terms appearing on the right band 
side of equation (lA.llI) when multiplied by Cr: 



' [ (c" + -C ] Crdrdt 
ti Jo \ r 



Ì2 



rC'C 



ti 



r=l 



dt 



r=0 



i2 fi 



ti JO 



{C'fdrdt 



(A.22) 



where again the boundary term vanishes. Finally, for the right band side of equation (IA.12I ) being 
multiplied by {£ - 2B) we use 



i2 /■! 



rdrdt 



t2 



-r£' (£ - 2B) 

k 



r=l 



r=0 



dt 



ti JO 



Ì2 



1 



Ì2 /"l 1 

-r {£') ^ - 2B'£'j drdt - j ^{£-2By 



r=l 



dt 



r=0 



(A.23) 



The first boundary term vanishes because £ satisfies Neumann conditions at r = 1 and manifestly 
vanishes at r = 0. The second boundary term does not vanish but since B satisfies Dirichlet conditions 
at r = 1 and the conical condition is satisfied at r = it equals 



Ì2 



-{£- 2Bf 

k 



r=l 



r=0 



dt 



Ì2 1 

-£'^ (t, l)dt<0 
, k 



(A.24) 



Note that without the conical condition imposed, the sign of this term would not be controlied and the 
uniqueness proof would fail. 

The other terms appearing on the right band side in the process of the multiplication and integration 
will be estimated as follows. For the difference of the exponentials we observe 



^4B+2C-2E 



AB+2C-2E 



dr L 



^{4B+2C-2E)t+{1-t)(4:B+2C-2È) 



dT 



^{4,B+2C-2E)T+{l-r)(iB+2C-2É) 

< c(|^| + |C| + |é:|) 







(4B + 2C- 2£) 



(A.25) 
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for some Constant C, following from the fact that the solutions {B, C, E) and {B, C, E) are bounded. 
Consequently, the lemaining teims on the right hand side of the B equation when multiplied by 
(j3r + Br^ can be estimated (note that the expression q'^b+2C~2E bounded above and below by 
some Constant) 



reniQ < K 



t2 fi 



ti JO 



+ + drdt + e 



t2 ri 



ti JO 



Here we frequently appUed Cauchy's inequahty 



ab < -a^ H . 

- 2 2i/ 



^^4B+2C~2E (^2\^^^^^ (A.26) 



(A.27) 



by means of which the e in (IA.26I ) can be made arbitrarily small on the cost of K becoming larger and 
larger. By the same token - since the remaining terms of the C equation are just given by differences 
of exponentials of the type considered in (IA.25I ) - we can estimate 



reme <K T f r {B^ + + S'^) 
Jti Jo 

For the remainder terms of the £ equation (multipUed by {£ — 2B)) we estimate 
1 



(A.28) 



k 



{£ - 2B) r 



C + C 



C + 2 



B + B 



B' 



2B + C-E + 2B + C-E 



{2B' + C' - £') 



<6-r U'^ + B'^ + C') +K-r{£^ + B^ + C^) 



(A.29) 



where we bave used inequality (IA.27I) and the fact that the spatial derivatives of the solutions are 
bounded. Taking the other terms (again differences of exponentials) into account we arrive at the 
estimate 



reme <K [ r {B'^ + + £^) drdt + 5 [ r (b'^ + C'^ + £'^) drdt (A.30) 
Jti Jo Jti Jo ^ ' 



t2 /•! 



where K is some Constant. Adding up the equations we fìnd the inequality 

r-l 



-r 



re 



AB+2C~2E 



B' 



dr 



t=t2 
t=ti 



+ - 



B'^rdr 



t=t2 



+ 



t=t^ 



t2 fi / 1 



4B+2C~2E ti2 



B'drdt 



< K 



t2 /■! 



il JO 



{B^+C^+£^) drdt{A3l) 
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We fix k first (say k = 16), then choose e, 5 small enough (say e = 5 = ^) to make ali terms on 
the left band side of expression ( IA.31I ) manifestly non-negative. This fixes the Constant K because 
(IA.27I ) has to be applied with a certain weight to produce the chosen e and 6 in the estimates (IA.26I ) 
and (IA.29I) . Since the time-dependent, manifestly non-negative expression 



P{t) = ^ I re 



1 

4:B+2C-2E 







dr 



1 



1 

2 



+ - / B' rdr (A.32) 

is zero at ti = we can find an interval [ti , i] with t < T in which expression (IA.32I ) is monotonically 
increasing or Constant in time. We then have 



max P{t)<K r [ riB^ + C"^ + drdt <K [ ^ P {t) dt < K {t2 - h) max P (t) 
^[tiM Jti Jo Jti te[ti,t2] 

(A.33) 



for some Constant K. Inequality ( IA.33I ) is valid for ali t2 < t. Finally, decomposing the time-interval 
[ti,t\ into a finite number of subintervals, each of size KAt = we find that P (t) = for each 
subinterval and since P (t) is manifestly non-negative, the solution has to be zero on each such subin- 
terval. But if i was not T, expression (IA.32I ) must be stiictly monotonically decreasing on a following 
small interval [i, t + e], which is impossible since it was just shown that the manifestly non-negative 
expression (IA.32I ) is zero alt = i. Hence t = T and B = C = £ = everywhere on [0, 1] x [0, T]. 

We note that a similar calculation can be performed for the equations following from the NUT-ansatz 

g = e2^(*''^)fir2 + r2e2^(*'-) {af + 4) + rh^^^'-'^al (A.34) 

with the gauge 

e = [- {2B' {t, r) + C (t, r) - A' (t, r)) + {2B'q (r) + (r) - A'^ (r))] dr (A.35) 

where AQ{r), Bq (r) and Co (r) are the functions chosen to obtain the Ricci-fiat Taub-NUT metric in 
(IA341) . 
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